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Abstract 

We study the propagation of sound in a homogeneous dipolar gas at zero temperature, known as 
zero sound. We find that undamped sound propagation is possible only in a range of solid angles 
around the direction of polarization of the dipoles. Above a critical dipole moment, we find an 
unstable mode, by which the gas collapses locally perpendicular to the dipoles' direction. 



Electronic address: bohii@murphy.colorado.edu 



I. INTRODUCTION 



At ultracold temperatures, the usual propagation of thermodynamic sound waves in a 
gas is suppressed due to the rarity of colhsions. Nevertheless, in a Fermi gas, collective 
excitations known as zero sound can still propagate even at zero temperature . These 
excitations involve a deformation of the Fermi surface in a way distinguished from that of 
ordinary sound. The theory of zero sound has been applied in the past to isotropic, two 
component Fermi liquids. Depending on the properties and strength of the interaction, one 
observes the possibility of different propagation modes (longitudinal or transverse) as well 
as the possibility of no zero sound propagation at all. Zero sound has been measured 
experimentally in liquid He-3 0, 4| in which both longitudinal and transverse modes of zero 
sound co-exist. Another type of collective excitations that exist for two component Fermi 
liquids is spin-dependent vibrations known as spin waves. 

Here we shall be interested in zero sound in a homogeneous, single component dipolar 
Fermi gas, where all the dipoles are oriented in one direction by an external field. In 
particular we are interested in the question of how zero sound propagation is affected by the 
long range and anisotropic nature of the dipolar interaction. In the case of Bose-Einstein 
condensates (EEC's), strong dipolar effects have been observed in a gas of Chromium-52 



atoms with relatively big magnetic moments 



,y,3- 



for a recent review of dipolar EEC's see 



. In the case of Fermi gases, the atomic magnetic dipolar interaction is small compared to 
the Fermi energy, so it may be hard to observe dipolar effects. However much stronger effects 
are expected in an ultracold gas of hetronuclear molecules with electric dipoles, where the 
dipolar interaction is considerably larger (such an ultracold gas of fermionic KRb molecules 
has been recently realized in experiment [o], 1^). Previous theoretical work on ultracold 



dipolar Fermi gases _s ;udied their ground state properties and expansion dynamics in the 



normal phase [11 



I2,ll3|,llj, as well as ECS superfluidity [l5|, ll6|. It has been shown that 



both the critical temperature for the superfluid transition, and the ECS order parameter, 
are sensitive to the aspect ratio of the trap. In the case of a cigar trap, the order parameter 
becomes non-monotonic as a function of distance from the center of the trap, and even 
switches sign. However, zero sound refers to phenomena in the normal phase, therefore we 
are interested in the regime of temperatures which are well below the Fermi temperature 
Tp, but still above the critical ECS temperature. Zero sound in Fermi dilute gases of alkali 



2 



type (with contact interactions), in the crossover from 3D to 2D, has been recently studied 



in Ref . 



17|. 



An interesting result due to Miyakawa eta al. ([12], uM) that the Fermi surface of 
a dipolar gas is deformed. They postulated an ellipsoidal variational ansatz for the Fermi 
surface, and by minimizing the total energy of the system, found that that it is deformed 
into a prolate ellipsoid in the direction of the polarization of the dipoles. As a preliminary 
to studying zero sound, one of our aims in the present work is to compute numerically 
the equilibrium Fermi surface, and thus, along the way, also check the reliability of the 
variational method. 

This paper is organized as follows. In section ([11]) we find numerically the shape of the 
Fermi surface of a dipolar Fermi gas at equilibrium, and compare it with the variational 
method. In section flIII|) we study the zero sound in this gas, and how it depends on the 
direction of propagation and the strength of the interaction. We also find the instability 
of the gas above a critical dipolar interaction strength, manifested in the appearance of a 
complex eigenfrequency. Finally we present our conclusions in section (1IVI) . 



II. SHAPE OF THE FERMI SURFACE 



We consider a homogeneous, single component, dipolar gas of fermions with mass m and 
magnetic or electric dipole moment d. The dipoles are assumed to be polarized along the 
2;-axis. It is significant that dipolar fermions continue to interact in the zero-temperature 
limit, even if they are in identical internal states. This is a consequence of the long-range 
dipolar interaction. The system is described by the Hamiltonian: 



1=1 i^j 



(1) 



where Vddir) = {d'^/r^){l — Sz'^/r'^) is the two-body dipolar interaction. We use the semi- 
classical approach in which the one-body density matrix is given by 12(] 



p(r, r') 



(2) 



{27ry \ 2 

where f{r,k) is the Wigner distribution function. The number density distributions in real 
and momentum space are given by 



n{r) = p{r,r) = J ^^f{r,k) (3) 
n(fc) = (2vr)"^y" d^rf{r,k). 

Within the Thomas Fermi-Dirac approximation, the total energy of the system is given 
by E = Ekin + Eii + Eex, the sum of kinetic, dipolar and exchange energies, where: 

. f d Ix, Ti }\, ^ . ^ . , . 



= i / <i'r / <i'r'14i(r - r')n(r)n(r'), (5) 



We shall consider a homogeneous system of volume V with number density nj. Let fcjT' 
be the Fermi wave number of an ideal Fermi gas with that density. Then is related io kp 
hj Uf = kp/^QTr"^). For an homogeneous gas the exchange energy can be rewritten as 

Here we have used the Fourier transform of the dipolar potential, Vdd{q) = 
{A7t/3)cP{3 cos^ 6q — 1) where 6q is the angle between the momentum q and the direction of 
polarization, which is chosen to be along the /c^-axis. 

For a homogeneous gas, the distribution function / is a function of k only. For an ideal gas 
or for a gas with isotropic interactions at zero temperature, it is given by f{k) = Q{kF — k), 
where Q{k) is Heaviside's step function. This describes a Fermi sphere with radius kp- 



In the case of dipolar gas, Miyakawa et al.([l4|, s[13|) postulated the following variational 

ansatz: 

m = eU-^{ki+ki)-p^kA, (8) 
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It is then possible to determine variationally the parameter (3 that minimizes the total energy 
of the system. In general it is found that /3 < 1, namely the Fermi surface is deformed into 
a prolate spheroid. This occurs due to the exchange interaction (Eq. ([7j)) being negative 
along the direction of polarization. 

We shall first be interested here in finding the accuracy of this variational method, and 
also derive some analytical results for the exchange energy, in a similar fashion to the well 
known exchange energy of a homogeneous electron gas. 

Instead of directly minimizing the total energy, our starting point is that the quasi-particle 
energy e[k) = ^^y-^^j-^ on the Fermi surface (i.e, the chemical potential) must be constant 
in equilibrium. The quasi-particle energy is given by: 

Further we shall assume that f{k) is either or 1, that is, there is a well defined and 
sharp Fermi surface. Our first calculation is perturbative. For small d, the solution of 
e(fc) =constant, is found. Let n be a unit direction vector in k-space. Let k = kn lie on the 
(deformed) Fermi surface. Then we have: 

k = n(^kF + ^J ^J{k')V,a{k,n - fc')) , (10) 

where, for weak interaction, / on the right side can be taken to be the distribution function 
of an ideal gas. The integral can be evaluated analytically by using the convolution theorem, 
with the result: 



m = kF + A!!^ (3 cos\e) - l) , (11) 

where Q is the angle between n and the 2;-axis. For small Eq. ffTTj) is consistent with a 
Fermi surface of slightly ellipsoidal shape, Eq. ([8]), with /? = 1 — A md^p ^ TMx^ confirms 
that the ellipsoidal ansatz is indeed correct for small d. The corresponding energy is then 

E/V 

and the chemical potential is 

5 



m l 207r2 4057r4^ M ' ^ ' 



h^kl / 28 md^kp 
2m I ~ 4057r2*-^ 



1-^(^)1- (13) 



Next, for arbitrary d, we solve numerically for e{k) =constant. The calculation is a self 
consistent iterative generalization of the perturbative calculation above, where the constant 
is adjusted to obtain a fixed number of particles. 

It is convenient to present our results in terms of a dimensionless parameter 

D ^ ^k,, (14) 

which signifies the strength of the dipolar interaction between two particles separated by 
a distance 1/kp, in units of the Fermi energy. This dimensionless expression appears in 
the equations for E/V and /i above. For reference, a gas density of 10^^ cm~^ with dipole 
moment of 1 Debye, and atomic mass of 100 amu, has D=5.8. In Fig. ([1]) we compare the 
Fermi surface found numerically with the variational Fermi surface for a strongly interacting 
dipolar gas with D = 12.0. We see that even for this relatively strong interaction, the 
variational method works quite well. The actual Fermi surface is not exactly an ellipsoid 
and is slightly less prolate. 

In Fig. ([2^) we compare the exact and variational results for the aspect ratio of the Fermi 
surface shape, and in Fig. ([2]d) we compare the chemical potentials. For both quantities the 
variational method generally works very well. 

Since our numerical results are close to those obtained by the variational method, we 
might also expect an instability, due to the inverse compressibility becoming negative, to 
occur around D ^ 12.5, as found in Ref. 1^. However, a cautionary note is in place, since 
for such strong interactions one can expect deviations from the mean field theory, which 
might stabilize the gas - as occurs in the case of a two-component Fermi gas with negative 
contact interactions 181 ]. 

Recently, it has been suggested jl^ that for strong enough dipolar interactions, bi-axial 
nematic phase should appear. Namely, above a critical interaction strength, the Fermi 
surface will distort into an ellipsoid with three different semi-axis, spontaneously breaking 
cylindrical symmetry. In our work, we assumed from the start that cylindrical symmetry 
of the equilibrium state is not broken. However, the bi-axial phase was predicted to occur 
around D ^ 30, which is above the compressional instability limit of D = 9.5. In what 
follows, we shall restrict ourselves to discuss zero sound for D < 9.5. In fact, we find below 
that the gas becomes unstable and collapses even below this limit. 
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FIG. 1: Shape of the Fermi surface of a homogeneous dipolar gas with interaction strength D = 
12.0. Sohd hue: numerical calculation. Dashed line: variational method 




D D 



FIG. 2: a) Aspect ratio of the Fermi surface of dipolar gas as a function of the interaction strength 
D. Solid line: numerical calculation. Dashed line: variational method, b) Chemical potential of 
dipolar gas as a function of the interaction strength D. fi is given in units of K^k'p/m. The 
numerical (solid) and variational (dashed) curves are virtually indistinguishable. 
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III. ZERO SOUND 



We now consider collective excitations of a homogeneous dipolar Fermi gas at zero temper- 
ature, i.e zero sound. The normal speed of thermodynamic sound (first sound) is = 
where P = — |^ is the pressure. However, thermodynamic sound is highly attenuated for a 
rarefied gas at low temperatures (in the coUisionless regime). The attenuation of sound is due 
to transfer of energy from the sound wave to random molecular motion (heat), and depends 
on the ratio of the collision rate to the sound frequency. Specifically, the coUisionless regime 
is obtained for ujt >> 1, where r is the mean time between collisions, and uj the frequency 
of sound. In particular, for a Fermi gas, the coUisionless regime is rapidly obtained when 
the temperature drops below the Fermi temperature, and collisions are quenched by Fermi 
statistics. However Landau [l| discovered that a collective excitation of a Fermi gas at zero 
temperature is still possible in this case. This zero sound is essentially a propagation of a 
deformation of the Fermi surface. The theory of zero sound as applied to a two component 
Fermi gas with short range interactions is described in standard textbooks, e.g P. fiojl. 

Our goal is to apply this theory to the case of the one component dipolar gas. The local 
quasi-particle energy at position r is given by 

eir k) = - (15) 

where the last term, the contribution of exchange interaction, has been calculated within 
local density approximation (as in Eq. ([9])), suitable for collective excitations with long 
wavelength. Note that, compared to Eq. ([9]), Eq. ( ITSl) shows an extra direct-interaction 
term involving Vdd{f — f')- For a spatially homogeneous distribution function this term 
vanishes, because the angular average of the dipolar interaction is zero. 
We apply a Boltzmann transport equation which reads: 



where /(/) is the collision integral, and e can be taken, to first order, to be the equilibrium 
quasi-particle energy. The coUisional integral can be neglected in the coUisionless regime we 
are interested in here. However, more generally, it is interesting to note that the dipolar 
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interaction is long range, giving rise to both a mean field interaction {Vdd{f — f') term in 
the right side of Eq. f|T5l) ) and a collisional cross section that enters /(/)• In this respect 
the dynamics of dipolar Fermi gas are an intermediate between, on the one hand, plasma 
dynamics which to a first approximation are controlled by the mean field potential (as 
in Vlasov equation), and on the other hand, the dynamics of Fermi liquids with short 
range interactions, which are controlled by the local exchange interaction and the collisional 
integral. 

For small deviations from equilibrium, the distribution function changes in the near 
vicinity of the equilibrium Fermi Surface, and can be written in the form 



f = fe + e^'^-^-'^v{k)5{e,{k) - /x), (17) 

where is the equilibrium chemical potential, and eg is as in Eq. ([9]). z^(fc) is a function 
defined on the equilibrium Fermi Surface, and will signify the eigenmode of the zero sound 
wave. This perturbation describes a mode with spatial wavenumber q and frequency u). 
Eqs. ( IT6!) and (ITS!) then give: 



(. - .(.) . ,) .(^) - -Hk) ■ . / dS'^^^§^^^ + (18) 



+vik)-qVMfdS'- 



u(k') 



{2TT)^h\v{k')\ 

where v(k) = In the above equation, k and k' lie on the equilibrium Fermi surface. 

Also the integrals are taken on this surface. The first integral is an exchange interaction 
term, and the second a direct term. They can be combined to obtain: 



r iVddiq) - Vdd{k - k')) u{k') 

The problem of zero sound of a dipolar gas is thus reduced to finding, for a given wavenum- 
ber q, the eigenvalues u and eigenmodes z/(fc). In fact, it can be seen from Eq. ( fT9|) that u 
is linear in \q\, and thus we expect a linear phonon spectrum with a speed of sound which in 
general depends on the relative direction of the propagation direction q with respect to the 
direction of polarization. It can be seen that for undamped vibrations, the speed of sound 
s{q) = ^ must exceed a critical velocity Wdamp(4) = niaxfc{t;(fc) ■ q}, where q = q/\q\. To 
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see this, replace v by another unknown function t'(fe) = (cu — v{\z) ■ q)v{h\ Then Eq. (fT9|) 
can be rewritten as: 



When Lo < maxfc{t;(A;) ■ q}, the integrand has a pole, which must be avoided by going 
around it in the complex plane, giving an imaginary part which signifies the decay of such 
a vibrational mode. 

Eq. (|T9l) can be solved numerically by discretizing it on the Fermi surface. We use an 
evenly spaced, 64x64 grid, in spherical coordinates 6, 0, where the 'north pole' is in the 
direction of polarization of the dipoles. In accordance with the discussion above, we select 
the undamped modes which satisfy s{q) > fdamp(9)- For that purpose, it is sufficient to look 
for a few eigenmodes with the largest eigenvalues, using an Arnoldi method. We checked 
our numerical code by calculating the modes for the case of an isotropic Fermi Fluid j^, 
and comparing to the known analytic solutions. 

Before discussing our results for the dipolar gas, it will be helpful to briefiy review the 
known results for two component isotropic Fermi liquids. The case relevant to us is that of 
same spin vibrations, in which the two components move in unison. For an isotropic fiuid 
the Fermi surface is a sphere. The interaction term 



which appears in Eq. (fT9|) . is replaced, in the case of isotropic liquids, by F{k,k') = 



ing. When Fi = and Fq > there is a longitudinal mode which is concentrated on the 
Fermi sphere around the forward direction of propagation, which becomes more and more 
concentrated around that point as Fq 0. For Fi = and Fq < 0, no zero sound propaga- 
tion exists. For Fi > 6, there also exists a transverse mode which has a vortex like structure 
on the Fermi surface around the direction of propagation. 

We now examine the results for the dipolar gas. Fig. shows the speed of sound s = u/q 
in the dipolar Fermi gas for interaction strength parameter D = 1, as a function of the angle 
a between the direction of sound propagation and the direction of polarization. We find 
one zero sound mode which, for small a, satisfies the criterion for propagation However 
sound does not propagate in all directions: for a > ~ 0.6 the gap between the calculated 




(20) 



Fiq,k,k') = VM-Vddik-k'), 



(21) 




where Fq and Fi are the Landau parameters. One finds the foUow- 
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speed of sound and the propagation limit Wdamp is essentially zero. This means that there 
is no undamped propagating zero sound for a wide range of angles around the direction 
perpendicular to the direction of polarization. 

To understand this result, we borrow lessons learned from the case of isotropic Fermi 
fluids. There, we know that spin-independent, longitudinal zero sound mode, only propa- 
gates for repulsive interaction (positive landau parameter Fq), and is damped for attractive 
interaction (negative landau parameter Fq). Thus, it makes sense that the anisotropy of 
the dipolar interaction gives rise to sound propagation only in a certain range of directions 
around the direction of polarization, where the effective interaction in momentum space is 
repulsive. 

For weak interactions, we can gain some insight into the observed behavior, as follows. In 
this case, the Fermi surface is nearly a sphere. Moreover, the eigenmode z^(fc) is appreciably 
different from zero only in a small region on the Fermi surface around the direction of q. This 
behavior is expected from the general Landau theory 2I and is confirmed by examination 
of our numerical solutions. Therefore, we may replace, in Eq. ( IT9l) . the interaction term 
F{q, k, k') with its value for k', k close to q. Since Vdd is anisotropic, F does not have in 
general a single-valued limit as k, k' q, because the limit depends on the direction of 
k' — k. Nevertheless we can observe the following. For the case a = 0, i.e, q along the 
direction of polarization, we find F does have the single valued limit of F = AtkP. This is 
completely analogous to a positive Landau parameter Fq in the isotropic Fermi liquid Q], 
and thus we obtain a propagating zero sound mode in this direction. On the other hand, 
for a = 7r/2 it is easy to estabhsh an upper bound F{q, k, k') < 0. Although here F does 
not have a single- valued limit, we can still expect, from analogy to the Fq < case in the 
isotropic Fermi liquid, that longitudinal zero mode does not propagate for this case. 

As observed above, numerically we determine the critical angle for propagation to be 
around a = 0.6 for D = 1. However we caution that we could not get to limit of very weak 
interactions D ^ hj our numerical method, since in this case the eigenmodes tend to be 
highly concentrated around the forward direction of propagation, and cannot be resolved 
with our grid. 

As we increase the interaction strength from D = 1 to D = 3 (Fig. (jl])), there is a 
noticeable change in the shape of the speed of sound curve. Yet, we still have one propagating 
mode with some critical angle beyond which zero sound does not propagate. 
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FIG. 3: Speed of sound (solid curve) in a dipolar Fermi gas as a function of the angle of propagation 
a relative to the direction of polarization, for dipolar interaction strength D = 1. The speed is 
measured in units vp = hkp/m. The dashed curve represents a lower bound on the speed of 
any undamped mode. When the two curves merge, zero sound modes are damped and fail to 
propagate. 

The deformation of the Fermi surface studied in the previous section contributes some- 
what to the change in the speed of sound. It somewhat changes the Fermi surface element 
and the Fermi velocity z^(fc) in Eq (fT9|) . However, the main factor giving rise to a change 
from propagation to damping of zero sound as a function of direction is the anisotropy of 
the dipolar interaction which directly enters that equation. 

For D > 4 we find that the solution of Eq.( [T^ admits complex eigenvalues, in particular 
an eigenvalue which is purely imaginary and lies on the upper complex plane. This signifies 
an instability, since according to Eq. (fT7|) such a mode grows exponentially. The imaginary 
eigenvalue appears first about D = 4.0 and a = 7t/2. Thus, the gas is unstable to growth of 
density waves in the plan perpendicular to the direction of the dipoles. For D > A we find 
also additional propagating modes (with real valued frequncies), in particular a transverse 
mode analog to the transverse mode in superfiuid He-3. However, because of the instability, 
we conclude that these modes could not be experimentally observed, at least unless some 
additional stabilizing mechanism is introduced. 

As noted in the previous section, the compressibility of the gas only becomes infinite at 
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FIG. 4: Speed of sound in a dipolar Fermi gas as a function of the angle of propagation a, for 
D = 3. 

D = 9.5. However it is physically reasonable that a dipolar gas becomes unstable prior 
to that. The reason is that the compressibility is a derivative of the pressure, which is 
an isotropic quantity. But due to the anisotropic interaction, the gas is actually more 
sensitive. to collapse in a specific direction, i.e, perpendicular to the dipoles. This tends 
toward the growth of over-density regions where the dipoles are oriented head to tail, lowering 
the potential energy. 

This mechanism is similar to that already long known for homogeneous dipolar Bose 
gases. A homogeneous Bose gas with a purely dipolar interaction is unstable, but it can be 
stabilized by the addition of repulsive short range interaction with scattering length a > 
{2^. In terms of a dimensionless dipolar interaction strength Db = md'^/h^a, the stability 
criterion of the homogeneous Bose gas is Db < 3. Also, from examining the dispersion 
relation for the Bose gas it can be shown that it becomes first unstable at = 3 due to 
density waves perpendicular to the direction of polarization. Finally, we note that Db and 
the dipolar Fermi parameter D (Eq. f[T^ are related by the simple exchange of the Fermi 
length scale l/kp with the scattering length a. It is interesting to observe that in these 
respective units, the criteria for the onset of instability in Bose and Fermi gases are very 
similar in magnitude. 

Experimental prospects: We outline some additional considerations regarding the exper- 
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imental prospects of observing zero sound in dipolar Fermi gas. The condition to have a 
normal phase is that T > Tc, where Tc, the critical temperature for the superfiuid transition, 
is given by ~ 1.44Tj exp (— 7r'^/(4Z})), with Tf the Fermi temperature [l^. At the same 
time, zero sound attenuation is proportional to T'^/t, where r is the relaxation time due 



to collisions [2l|]. This will set some maximal temperature for realistic experimental 
detection of zero sound. Determining the zero sound attenuation at finite temperature the- 
oretically is a task beyond the scope of this paper , but we note the following: observing 
zero sound in the normal phase requires > T > T^, and we thus obtain the condition 
D < 7rV(41og(1.44T//r„)). As an illustration, if = O.IT/, we obtain D < 2.9. Still, 
even for smaller T^, restricting the observation of the zero sound to smaller D, it should 
be possible to see the propagation limited to certain directions only (as in Fig. ([3]). Finally, 
the above discussion is based on the requirement that T > Tc. Nevertheless, we note that, 
under certain conditions, zero sound of the normal phase can still persist below T^, into the 



superfiuid phase |22l. |23|. 



IV. CONCLUSIONS 

In conclusion, we have studied numerically the deformation of the Fermi surface and zero 
sound excitations in a homogeneous, single component, degenerate Fermi gas of polar atoms 
or molecules polarized by an external field. We find that the Fermi surface is described very 



well by the variational ansatz proposed in 12(] up to about dipolar interactions strength 



D ^ 12.5, where a compressional instability was predicted. We then study zero sound in 
the range D = — 10. For D < 4 we find that zero sound can propagate with a longitudinal 
mode in directions between a = and some critical a which depends weakly on D but is 
generally close to 0.6 — 0.7 radians {a being the angle between the direction of propagation 
and the direction of polarization). Beyond this angle, undamped propagation of zero sound 
is not possible. For D > 4, we find a complex eigenfrequency indicating that the dipolar gas 
collapses due to an unstable mode in the direction perpendicular to the dipoles. 



Note added: after completion of this work, a work by Chan et al. appeared 
study properties of dipolar Fermi gas, including zero sound. 



24| . who also 
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